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ABSTRACT 

We present high spatial resolution Chandra X-ray images of the NGC 2237 
young stellar cluster on the periphery of the Rosette Nebula. We detect 168 X- 
ray sources, 80% of which have stellar counterparts in USNO, 2MASS, and deep 
FLAMINGOS images. These constitute the first census of the cluster members 
with 0.2 < M < 2 Mq. Star locations in near-infrared color-magnitude diagrams 
indicate a cluster age around 2 Myr with a visual extinction of 1 < Ay < 3 at 
1.4 kpc, the distance of the Rosette Nebula's main cluster NGC 2244. We derive 
the K-band luminosity function and the X-ray luminosity function of the cluster, 
which indicate a population ~400-600 stars. The X-ray-selected sample shows a 
i^-excess disk frequency of 13%. The young Class II counterparts are aligned in 
an arc ~ 3 pc long suggestive of a triggered formation process induced by the O 
stars in NGC 2244. The diskless Class III sources are more dispersed. Several 
X-ray emitting stars are located inside the molecular cloud and around gaseous 
pillars projecting from the cloud. These stars, together with a previously unre- 
ported optical outflow originating inside the cloud, indicate that star formation 
is continuing at a low level and the cluster is still growing. 

This X-ray view of young stars on the western side of the Rosette Nebula 
complements our earlier studies of the central cluster NGC 2244 and the em- 
bedded clusters on the eastern side of the Nebula. The large scale distribution 
of the clusters and molecular material is consistent with a scenario in which the 
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rich central NGC 2244 cluster formed first, and its expanding HIT region trig- 
gered the formation of the now-unobscured satellite clusters RMC XA and NGC 
2237. A large swept-up shell material around the HII region is now in a second 
phase of coUect-and-collapse fragmentation, leading to the recent formation of 
subclusters. Other clusters deeper in the molecular cloud appear unaffected by 
the Rosette Nebula expansion. 

Subject headings: Open clusters and associations: individual (NGC 2237, NGC 
2244) - ISM: individual (Rosette Nebula) - stars: formation - stars: pre-main 
sequence - X-Rays: stars 



1. Introduction 



The triggered formation of the lower mass stars in the vicini ty of massive stars is a 
complex process that is only now being characterized in detail (see iBriceno et al.l 120071 for 
a recent review). In their immediate neighborhood, massive stars suppress f urther star for - 
mation by quickly ionizing and dispersing surrounding molecular material (iHerbigl Il962l ). 
At greater distances, OB stars are more constructive to star formation activity; the shocks 
driven by ionization or stellar winds are crucial i n triggering the collapse of molecular cores 
( iWhitworth et al.lll994l : iLefloch fc Lazareflall994j ). Triggered star formation events by mas- 
sive stars have been observed at diffe rent spatial scales, for example, sraall bright-rimmed 



clouds on the periphery of HII regions ( ISugitani et al.lll995l : iGetman et al.ll2007l : lOgura et al. 



20071). an embedded cluster in a molecular cloud core on the edge of the dispersed Cep 0B3 b 



( iGetman et al.l 120061 ). multiple generations of star formation in W5 (iKoenig et al. 



a broad ridge of young stars along the southwestern boundary of M 17 ( 



Broos et al. 



20061 ). 



20071 ). and a rich secondary cluster on the edge of Sharpless 219 (jPeharveng et al. 



iang et al. 



20081), 



20021: 



The Rosette star forming complex has been considered an excellent candidate for trig- 
gered st ar formation (ICox et al.lll990l : iPhelps fc Ladal 119971 ) following the framework devel- 
oped W Elme^^een^&Lada (jl977l ). The massive young cluster NGC 2244 {d ~ 1.4 kpc, t ~ 2 
Myr; iHensberge et al.ll2000l ) powers a visually spectacular expanding HII region known as the 
Rosette Nebula. The ionized nebula is clearly interacting with the adjacent Rosette Molecu- 
lar Cloud (RMC) to the east of the Nebula which has a collection of embedded young stellar 
clusters, each with a few hundred pre-main sequence stars (Table 6 in Wang et al. 2009; see 
also Poulton et al. 2008). However, recent investigations of the ages and disk fractions of 
these embedded stellar populations at near- infrared (NIR), mi d-infrared, and X-ray wave- 
lengths do not obviously support a sequential, triggered origin ( iRoman-Zufiiga et al.ll2008al : 
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Poulton et al.l l2008l : IWang et al.l 120081 ) . Instead, Wang et al. suggest that the previously 
unnoticed, more evolved, and less-obscured RMC XA cluster lying between NGC 2244 and 
the molecular cloud may have been triggered by the expanding Rosette HIT region ~ 2 Myr 
in the past. 

To the west of NGC 2244 cluster hes the NGC 2237 clusteiQ. NGC 2237 was not studied 



until the recent 2MASS and FLAMINGOS JHK surveys ^ l2005l : iRoman-Zuniga et al. 
2008al : iBonatto fc Bicall2009l ). Li (2005) estimates th at the cluster has a diameter of 3.2 pc 
with a population of 232 stars. ISonatto fc Bkal J2009h studied NGC 2237 (called NGC 2239) 
with field-star-decontaminated 2MASS photometry, and suggested that it may be a young 
(5 ± 4 Myr) cluster located in the background of NGC 2244 {d = 3.9 kpc). In the Digital 
Sky Survey optical plate (Figure [T]), the NGC 2237 region lies o n the boundary between 



the HII region and the cold n eutral molecular material (Core E; iBlitz &: ThaddeusI Il980 



Roman-Zuiiiga fc Ladall2008bl ) on the western side of the nebula. Numerous dusty elephant 
trunks and pillars point towards the O stars in NGC 2244. Its structure contributes to the 
petal-like morphology of the Rosette Nebula and resembles the pillared boundaries of other 
HII regions such as M 16 and IC 1396 where the molecular cloud is photoionized and ablated 
by the O star ultraviolet light and winds. 

However, the stellar content of NGC 2237 has not been well-characterized and most of 
its members have not been individually identified. The difficulty is that the near-infrared 
surveys are dominated by older field stars unrelated to the Rosette star forming complex so 
that only the statistical enhancement in stellar surface density, and the fraction of members 
with i^T-band excesses, have been studied. X-ray surveys are more effective at uncovering 
the full pre-main sequence (PMS) population of young clusters. Low mass young stars 
emit X-rays from violent magnetic reconnection fl ares, orders of magni tude more luminous 
than seen in older Galactic stars (see review by iFeigelson et al.l 120071 ). Hence, an X-ray 
study can individually identify young stars including the dis kless Class HI sou rces missed 
in the infrared-excess samples. X-ray studies of NGC 2244 (jWang et al.ll2008l henceforth 
Paper I) and the embedded clusters in the RMC (IWang et al.ll2009l . Paper II) on the eastern 



^The designations of clusters and nebulae in the Rosette region have been confused since the 19th century. 
NGC 2 237 refers both to the stellar cluster and its associated nebulosity on the west side of the Rosette 
nebula ( Sulentic fc Tifft 1973 ). The star cluster is listed as O cl 512 i n the catalogue of Alter et al. ( 19701 ) 
and included in the study of the Monoceros star clusters by iPerea (jl99l[ ). NGC 2238 and 2246 are old 
designations for bright emission regions around the Rosette Nebula; these names are rarely used today. The 
well-studied rich star cluster exciting the nebula is NGC 2244; NGC 2239 is an obsolete designation for this 
central cluster. The listing in the SI MBAD database f or NG C 2237 was incorrect until recently. Note that 
the open cluster NGC 2239 studied in lBonatto fc Bica (l2009h in fact refers to NGC 22 37. A thorough review 



of the Rosette star formation complex is provided by 



Roman- Zuhiga &: Lada (l2008bl) 
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side of the Rosette Nebula produced a richer census of young stars than available from 
infrared-excess studies. A handful of X-ray stars in NGC 2237 have been identified with 



the ROSAT satellite (see iBerghofer fc ChristianI l2002l ) . but observations with the Chandra 



X-ray Observatory with its high spatial resolution and sensitivity, are much more effective. 

In this paper (the third in a series on the Rosette Complex), we focus on the X-ray point 
sources identified in the image of the NGC 2237 region obtained with the Chandra X-ray 
Observatory. We identify over 150 individual X-ray emitting young stars and investigate 
their spatial distribution and stellar properties to study the star formation activity on the 
western side of the Rosette Nebula. Proximity to dense molecular structures and correlated 
age-spatial gradients in the stellar distributions are particularly valuable in evaluating the 
timescale and efficiency of the triggering process that may have occurred. We combine these 
results on NGC 2237 with earlier work on the central NGC 2244 cluster and the embedded 
RMC clusters to elucidate a large-scale view of the star formation history in the Rosette 
complex. A distance of d = 1.4 kpc to NGC 2237 is adopted throughout the paper and the 
uncertainty in distance (see lBonatto fc Bical 120091 ) will be discussed. 



2. Chandra Observations and Data Reduction 



The NGC 2237 region was observed with the Imaging Array of the Chandra Advanced 
CCD Imaging Spectrometer (ACIS-I) which has a field of view ~ 17'xl7'. The 20 ks ACIS-I 
observation was conducted on 9 February 2007, in the "Timed Event, Very Faint" mode 
with 5 pixel x 5 pixel event islands (see Table 1 in Paper I for details of pointing and roll 
angle). The pointings towards NGC 2237 and NGC 2244 are outlined in Figured! 

We follow the data reduction and source extraction procedures described in Paper I. 
The ACIS-I image is shown in Figu re [21 A smoothed, exposure corrected X-ray image is 
created with the CIAO tool csmooth (lEbeling et al.ll2006l ) and shown in Figure [3l which gives 
the relative spectral hardness of the sources. 65% of the sources are dominated by soft band 
(0.5-2.0 keV) X-ray emission with median photon energy <2 keV, implying re latively low ob- 



scura tion. A candidate source list was obtained using the wavdetect program ( iFreeman et al. 



2OO2I ). supplemented a few possible sources found by visual inspection. Events are extracted 
and analyzed for each candidate source using our IDL script ACIS Extract ( version 3.128; 
hereafter AE, iBroos et al.ll2010l ). After a careful review of the net counts distribution for 
all candidate sources, we rejected sources with greater than 1% likelihood of arising from a 
background fiuctuation assuming Poisson statistics {Pb > 0.01). Only the brightest source 



^http : //www . astro . psu . edu/xray/ docs/TARA/ ae_users_guide . html 
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in the field (#149) has more than 100 source counts. 

The trimmed source hst includes 168 valid X-ray detections, which are further divided 
into a primary list of 130 highly reliable sources (Table [T]) and a secondary list of 38 tentative 
sources with 0.001 < Pb < 0.01 likelihood of being spurious background fluctuations. The 
resulting sources and X-ray properties are reported in Table [J a nd [2] which are very similar 
to the source tables in lTownsley et al.l (120061 ) . iBroos et al.l (120071 ) . and Paper I. Descriptions 
of the table columns are given in the table footnotes. Six sources (#160, 163, 164, 166, 167, 
and 168) are in the ACIS region overlapping the NGC 2244 field described in Paper I. Four 
of these (#163, 166, 167, and 168) were detected in the earlier ACIS ex p osure . One bright 
source (#149) was found in the i?(95'^r study of iBerghofer &: ChristianI (120021 ). 



Source event arrival times were compared to that of a uniform light curve model in 
AE, and five sources (#4, 39, 118, 149, and 152) were found to be significantly variable 
{Pks < 0.005 in column 15 of Tables [1] and [2]) during the short 2 ks observation. The light 
curve of #149 shows a typical X-ray flare seen in PMS stars (e.g. IWolk et al.ll2005l ) captured 
in the rise phase near the end of the observation. 

For the 65 sources with photometric significance Signif > 2.0 (column 12 in Tables [1] 
and El) , the extract ed spectra were fit with a single temperature apec thermal plasma model 



( Smith et al.ll200ll ) subjected to an absorbing column of interstel lar material ([Wilms et al. 



2000|, tbabs model; ) with the X5PE(f| package (version 12.2. lap lArnaudlll996h . The best- 
fit model was obtained with the maximum likelihood method (jCashl Il979l ). Abundances 
of 0.3 Zq were assumed for the automated fitting performed by AE. Spectral results are 
presented in Table [3l table notes describe the columns and unusual cases. Absorbing column 
densities range f rom negligible to lo gA^ff ~ 23.2 cm~^, equivalent to a visual absorption of 
Ay ~ 100 mag (jVuong et al.ll2003l ). Temperatures range from kT ~ 0.1 keV to kT > 15 
keV; this truncated value reflects the hardest spectra detectable by the ACIS detector. The 
range of total band (0.5 — 8 keV) absorption-corrected luminosities derived from spectral 
modeling is 30.0 < log L^^c ^ 32.1 ergs s~^. 

To estimate the X-ray sensitivity of the ACIS exposure, we use the Portable, Interac- 
tive Multi-Mission Simulator (PIMMS) softwar^ for a limiting on-axis source with 3 counts 
assuming a 2 keV plasma temperature and an absorbing column log Nh ~ 21.4 cm~^ (corre- 
sponding to Av = 1.5 mag visual extinction). The resulting limiting total band luminosity 
as observed by the telescope is logL^ ~ 29.4 ergs s~^ and limiting intrinsic luminosity 



http: / /heasarc. gsfc.nasa.gov/docs / software /Iheasoft / xanadu /xspec 



http : //heasarc . gsf c . nasa . gov/ docs/ software/ tools/pimms . html 
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log Lt^c ~ 29.6 ergs s~^. A conservative estimate for limiting luminosity of the entire NGC 
2237 observation, including sensitivity degradation due to off-axis optics and higher absorp- 
tion, is log Lt ~ 30.0 erg s~^. Using the empirical correlation between X-ray luminosity and 



mass in PMS stars from we ll-studied samples in the Orion and Taurus clouds (jPreibisch et al. 



20051 : iTelleschi et al.ll2007l ). we infer that the observation is complete to all cluster members 
down to M ~ 1.0 Mq and includes an incomplete fraction of cluster members down to 
M ~ 0.2 Mq. 



3. Stellar Counterparts to Chandra Sources 

3.1. Correlation of X-ray and Optical/Infrared Catalogs 

We searched optical and infrared catalogs from the literature and recent observations for 
stellar counterparts to these 168 X-ray sources. In addition to the available catalogs listed 
in Paper II (USNO, 2MASS, and FLAMINGOS), other catalogs are: UBV photometry of 
NGC 2244 (Ogura & Ishida 1981, 0181; Massey et al. 1995, MJD95), BVIRB.^ photometry 
of NGC 2244 (Berghofer & Christian 2002, BC02). Spitzer Space Telescope mid-infrared 
observations are not yet available for this field. The reference frame offsets between the 
ACIS field (which is astrometrically aligned to the Hipparcos frame using 2MASS sources in 
the data reduction) and the catalogs are 0.3" to 0181, 0.4" to MJD95, 0.3" to BC02, 0.2" to 
USNO, and 0.2" to FLAMINGOS. These offsets were applied before matching sources. 

We associate ACIS X-ray sources with optical and near IR (ONIR) sources using po- 
sitional c oincidence based on the positional uncertainties of the X-ray sources and stellar 



catalogs (IBroos et al.l 120071 ). Likely associations for 134 of the 168 ACIS sources are given 
in Table m along with ONIR photometric data (see Table H] footnote for photometric errors). 
Typical offsets between the X-ray and near-infrared sources are 0.2" for on-axis sources and 
0.7" for off-axis sources. Although USNO photometry is reported when optical counterparts 
are found, readers are cautioned that m agnitudes are deriy ed from photographic plates with 



0.3 magnitude photometric accuracy ( iMonet et al.ll2003l ). JHK magnitudes from 2MASS 



are given when the source falls outside the FLAMINGOS field of view. The SIMBAD and 
VizieR online catalogs are searched for complementary information with results given in the 
table footnotes. 

The Chandra observati on is a ble to identify a majority of the estimated ~ 230 members 



of the NGC 2237 cluster (jLill2005l ). This excellent result arises from a combination of three 
effects: the short ACIS exposure is sufficient to detect a large fraction of the stellar initial 
mass function (IMF), even with significant obscuration (^; the catalog from FLAMINGOS 
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deep JHK imaging is sensitive to most of the cluster population; and the X-ray image suffers 
only minor contamination from non-cluster X-ray sources (§ 13. 3p . 



3.2. Near-infrared and Disk Properties 



The J — H vs. H — K color-color diagram for 119 Chandra stars is shown in Figure |H 
This sample is restricted to stars with high-quality JHK photometry (i.e., errors in both 
J — H and H — K colors < 0.1 mag). Most Chandra sources occupy the color space associated 
with diskless young stars, Class III or weak-lined T Tauri stars (WTTS), that are reddened 
by interstellar extinction. A concentration of cluster members subjected to Ay ~ 2 mag is 
apparent, approximately centered at J — H = 0.80 and H — K = 0.35. This is also the typica l 



extinction seen towards NGC 2244 (e.g., lOgura fc Ishidalll98ll : lBerghofer fc Christianll2002l ). 
Seven stars with i^-band excesses lie to the right of this reddening bancjfl, defined as stars 
that have colors {J - H) > 1.7{H - K) + 2a{H - K). These are likely PMS stars with 
circumstellar accretion disks. Class II objects or classical T Tauri stars (CTTS). One outlier 
among the Class III sources, #86, appears highly obscured with Ay ~ 20 mag. Among the 
faintest X-ray sources with only 4 net counts in the ACIS image, this source is located in 
an optically dark region ; this so urce is discussed in § 4.2. We did not detect the red sources 



m 



Roman-Zuniga et al.l (l2008al ) with no detection in J and H — K > 1.5 mag, the majority 



of which were located within the molecular cloud (cf. Figure 27 in Roman-Ziifiiga et al.). 
However, four ACIS sources (#45, 70, 73, 85) have only i^-band counterparts with no J 
and H detections, and hard median photon energy (~2.5 keV). They are probably deeply 
embedded young stars, similar to those red sources. 

Figure shows the NIR J vs. J — H color-magnitude diagram (CMD) for the same 
stars pre sented in F i gure HI The 1 Myr and 2 Myr PMS isochrones from the evolutionary 
tracks of ISiess et al.l (120001 ) are shown for reference. If the NGC 2237 satellite cluster was 
triggered by the centra l NGC 2244 cluster, it must be younger than NGC 2244 at ~ 2 Myr 
(IHensberge et al.ll2000l ). While it is not possible to directly measure the age of NGC 2237 
without obtaining spectra of the cluster members, the distribution of stars in NGC 2237 in 
the NIR color-magnitude diagram is nearly identical to that seen in the main NGC 2244 
cluster (cf. Figure 6 of Paper I). Approximate stellar masses can be measured from the iso- 
mass vectors (blue dashed lines) on the color-magnitude diagram. They range from early-B 



■^These sources are #44, 61, 80, 97, 111, 117, and 121. For the two sources (#111 and #121) that are 
close to the reddening band, the errors in both J — H and H — K colors are small, ^ 0.015 mag, therefore 
their ii'-excesses are small but statistically signi ficant. See also Paper I for a brief discussion on systematics 
in the ii'-excess. These stars are also present in lRoman-Zufiiga et al.l (|2008al) . 
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stars to early-M stars. 

From Figures m and we find that the fraction of K-band excess sources is 13% among 
stars with M > 1 Mq, where we estimate masses from the NIR color-magnitude diagram 
(recall this is also the completeness limit of the X-ray-selected stars). This closely matches 
the K-excess disk fraction of 15 ± 3% reported in the FLAMINGOS NIR study of NGC 



2237 (IRoman-Zuniga et al.l l2008al ). However, when the full X-ray sample including many 
sub-solar mass stars is considered, the K-hand excess fraction of the X-ray sample is 6%. 
Since disks around intermediate-mass stars evolve faster than disks around low-mass stars 
(e.g., Hernandez et al. 2007), most likely this is a selection effect against the X-ray emission 
from CTTS which are known to be systematically lower than the X-rays from WTTS by a 
factor of 2 (e.g., Preibisch et al. 2005). 



3.3. Non-cluster Contaminants and Embedded Sources 



We deduce from previous studies (e.g., Wang et al. 2007) that the 34 ACIS sources with- 
out any matched counterparts are a mixture of newly discovered embedded cluster members, 
distant background stars, and extragalactic sources. The level of contamination by extra- 
galactic X-ray sources and Galactic disk stars is evaluated with simulations as described 
in Paper I. The extragalactic population is simulated by placing artificial sources randomly 
across t he detector with inc ident fluxes drawn from the X-ray background log A^-log S distri- 
bution (IMoretti et al.ll2003l ). The Galactic stellar population expected ir i the directio n of ou r 
ACIS field is simulated with the stellar population synthesis model of iRobin et al.l ( l2003l ). 
adopting X-ray lurn inosity distributions of nearby stars measured from ROSAT surveys (e.g., 
Hiinsch et al.lll999l ). These simulations predict that ~ 12 sources will be extragalactic and 
~10 could be old field stars. These contaminants constitute ~ 13% of the 168 ACIS sources. 
Figure [7] compares the spatial distribution of X-ray sources without NIR counterparts and 
sources that have hard median photon energies {Emed > 2.0 keV). The 10 unidentified X-ray 
sources that are located in the unobscured HII region cavity should have counterparts if they 
are cluster members. These sources, plus a few in the obscured region, are most likely the 
background extragalactic sources and background field stars. The other 24 ACIS sources 
lie in the molecular pillars and the optically darker regions with dense molecular material. 
They are likely new embedded low-mass PMS stars that are fainter than the FLAMINGOS 
completeness limit. Such X -ray-discovered stars are commonly foun d in the molecular clouds 
surrounding other clusters (ICetman et al.ll2005al : iBroos et al.l 120071 . Paper I). 
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3.4. Properties of Selected Stars 



A few stars have bright J magnitudes indicating they are intermediate-mass stars (B 
and A types) reddened from the ZAMS with < Ay < 2 mag. Their locations in the color- 
color diagram (Figure S]) are also consistent with early spectral types, and their bright optical 
counterparts are identified by 0181, MJD95 or other studies. Two of these sources (#44 and 
#61) have masses M > 2Mq estimated from the color-magnitude diagram and strong K- 
band excesses; these are good candidates for being accreting Herbig Ae/Be (HAeBe) stars. 
Others do not show NIR excess. Wang et al. (2007) also find that in massive OB cluster 
NGC 6357, over 90% of the X-ray selected intermediate-mass stars a re not HAeBe stars. 
Perh aps this is related to the fast disk evolution in the HAeBe stars (e.g.. lWaters fc Waelkens 
19981) . or an X-ray defic iency in detecting the intermediate stars with NIR-excess disks (e.g., 
Hamaguchi et al.ll2005[ ). 



ACIS source #53 is a Bl star (jPhilip fc Egretlll980[ ). the earliest spectral type known 
in the NGC 2237 region. It only h as 16 net counts and t he spectral fitting gives logLj c ~ 
30.2 erg s~^. A B2 star identified by lPhilip fc Egretl (Il980l ) is source #54 in our X-ray study. 
With over 100 counts, its X-ray spectrum can be reasonably well-determined (Table [3] and 
Figure |6^): it is dominated by a plasma with moderate temperature {kT = 1.2 keV) with 
low absorption (logNH ~ 20.0 cm~^) and X-ray luminosity logL^ c ~ 31.0 erg s~^ (0.5-8 
keV). These X-ray pro perties are similar to those of other early-B stars; for example, in the 
Orion Nebula Cluster dstelzer et all 120051 ONC;) and NGC 2244 Rosette Cluster (Paper 1). 



Source #149 is the X-ray source with most net counts in this study with a spectrum 
well- fit by a single plasma of kT = 4.1 keV and absorption column of logA^j^ ^ 21.3 cm~^ 
(Figu re [Hb)- The kT value is higher than that of a typical PMS star (e.g., iGiidel et al. 
20071 ). caused by the X-ray fiare seen near the end of the observation. It is known from pre- 
vious obser vations that spectral hardening is commonly seen during X-ray fiaring in PMS 



stars (e.g., iGetman et al.ll2008al ). The absorption corrected tot al band (0.5-8 keV) X-ray 



luminosity logLj^ ~ 31.3 erg s ^. iBerghofer fc ChristianI ( 12002| ) reported a similar level of 
logLj; = 31.01 erg s~^ (0.5-2 keV) in their 05*^ T/PSPC observation. This source repre- 



sents o ne of the typical fiaring PMS stars se en in the ONC (jWolk et al.ll2005l : iGetman et al. 
2008al ) and many other young star clusters (IWolk et al.ll2008l ). 



Lastly, ACIS source #20 is very interesting in its X-ray properties, showing a high X-ray 
luminosity (logL^c ~ 31.3 erg s~^), hard spectrum [kT ~ 3.4 keV), and high absorption 
{Nh ~ 10^^ cm~^). It is located in the optically dark pillar region at the peripheries of the 
HII nebula. No counterpart is found in the 2MASS image and the FLAMINGOS image, and 
no information on its mid-IR or mm band properties is currently available. It may represent 
a similar detection of Class O/I protostar as the VLA/BIMA source IRAS 21391-^5802 in IC 
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1396N (Getman et al. 2007) or a chance superposition of a background AGN. 



4. Cluster Properties 



4.1. Cluster Center 



The centr al posi tion of the NGC 2237 cluster has not been well- determined from NIR ob- 



servations. 



LJ (120051 ) obtain s (a,6) = (6^30"^55■8^ +04°5 8^30") from a weak 2MASS surface 
density enhancement, while iRoman-Zuniga et al.l ( l2008a[ ) obtains (6^30™' 22^, -t-04°560 from 



the -ft^- band excess sources in the deeper FLAMINGOS data (see Table 1 in iRoman-Zuniga et al, 
2008al ). The 2MASS position is ~8' east of the FLAMINGOS position. 



The X-ray position for the cluster core and the spatial distributions of stellar populations 
are best examined from a smoothed map of the X-ray selected cluster members. Following 
our Papers I and II, we derived a smoothed map of the surface density of the X-ray-selected 
cluster members using a 1.5'-radius top-hat smoothing kernel (Figure [8]). To reduce the 
effect of sensitivity variations across the field, we included only 93 sources that have 6 or 
more net counts. Adopting different kernel widths does not have a significant effect on the 
resulting features discussed here. The stellar density in the western portion of the field 
with molecular cloud material is likely underestimated compared to the eastern portion for 
two reasons: 24 unidentified X-ray sources are probably embedded stars (§ 13. 3p and the 
sensitivity to embedded stars is reduced due to soft X-ray absorption. 

A clear peak stellar concentration appears at (6'*31™00'^, +04°58'). This is the most ac- 
curate available as the X-ray contour map has a much higher dynamic range due to its smaller 
contamination by extraneous sources. This position derived from the Chandra selected stars 
i s con sistent with the low-resolution ROSAT X.-YaY image shown by iBerghofer fc Christian 
(120021 ). and agrees reasonably well with 2MASS result, lying 1' towards the main NGC 2244 
cluster. The discrepancy in the infrared cluster centroids is due to different fields of view, 
definition of cluster center, and possibly age effects: t h e FLA MINGOS field extends ~ 10' 
westward of the ACIS field, and iRoman-Zuniga et al.l (l2008al ) choose to define the cluster 
to encompass Class II systems found in this extended region, while the X-ray map and the 
2MASS K surface density are weighted toward Class III sources. Our X-ray star distribution 
shows a secondary peak 8' northw est of the primary peak (see 14. 2p . which is close to the 
NGC 2237 core position defined in iRoman-Zuniga et al.l (l2008al ). 
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4.2. Spatial Distribution of the Young Stars 

The spatial distribution of cluster members has not been well-established from NIR stud- 
ies due to the high contamination level of non-cluster stars. The diameter of the NGC 2237 



cluster has been estimated to be ~ 3 .2 pc (8') from the 2MASS e nhancement (lLill2005l ) and 



3.8 pc from the i^'-band excess stars (IRoman-Zuniga et al.ll2008al ). Our X-ray-selected NGC 



2237 stars, which represent the bulk of cluster members, indicate that the structure is more 
complex than a spherical aggregate with a single size. 

From Figure [H we see that the dimension of the main stellar concentration around the 
cluster core at (6'^31™'00*, +04°58') is ~ 6' x 3' (~ 2.5 x 1.3 pc) in diameter, asymmetrically 
elongated towards the west from the dense core. A sparse secondary stellar concentration 
may also be present: a northwestern group of ~ 15 stars in a 1.5 x 1 pc region around 
(6'*30™39^ +05°00'). 

We next investigate whether any spatial segregation is seen in stars at different stages 
of PMS evolution. Figure [9] shows the distribution of X-ray-selected Class II and Class III 
sources classified using the near- infrared color-color diagram (Figure H]). The 34 sources 
without matched counterparts, many of which are probably embedded stars (§ 13. 3p . are also 
shown. The seven Class II sources appear to be aligned along an arc parallel to the edge of 
the dusty cloud edge and perpendicular to the direction towards the early O stars in NGC 
2244. Most of the Class III sources are associated with the groupings mentioned above; the 
southern and eastern edges of the field seem to lack X-ray sources. The CO emission contours 
suggest that there is some molecular material to the west. The sensitivity of our observation 
is too shallow to probe the star formation activity here because of heavy absorption. 

Figure [10] shows the spatial distribution of stars in different mass ranges. We consider 
stars with identified B and A spectral types (see Table H] footnote) and stars with NIR- 
estimated masses > 2Mq as the high mass sample, while stars in the mass range M < 2Mq 
are the low mass sample. The high mass stars are spatially concentrated at the cluster 
center. Note that sources #44, 53, and 61 lie along the same arc as the i^-band excess Class 
II sources. In contrast, the majority of the low mass stars are not strongly clustered, but 
are dispersed along the peripheries of the HII region and inside the dark cloud. Again, the 
deficit of sources in the southern and western edges of the field is worth noting. 



A close-up Ha image (iLi &: Rectorl 120041 ) of a small region at the interface between the 



HII region and the cloud is shown in Figure [TTl A number of X-ray stars are found at or near 
the tip of dusty pillars (ACIS #77, 86, 91, 96, 100, and 101) in the two regions outlined by 
circles. The Chandra images of the Eagle Nebula (M 16) and NGC 6357 are similar with a 
handful of hard-spectrum X-ray stars embedded in dusty pillars on the edges of HII regions 
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ionized by rich clusters (Linsky et al. 2007, Wang et al. 2007). It is plausible that these stars 
become exposed because the OB stars of NGC 2244 (to the east) have photoevaporated the 
surrounding cloud material. The most obscured source, #86, is embedded in or behi nd the 



southern pillar shown. Only 4 counts are seen from this source, most above 2 keV. iLj (120051 ) 
also noted a highly obscured 2MASS source with extinction Ay ~ 17 mag located behind 
the same pillar object, however no matching X-ray source was detected in our observation, 
likely due to the limited sensitivity. 

Another interesting feature in this region is a bright visually identified outflow outlined 
by the white box in Figure [TTl Located at (6'^30'^48.7'*, +05°01'44"), it has not been reported 
previously in the literature, to our knowledge. This may be a Herbig-Haro (HH) object. 
There is a visible star 30"/ to the northwest aligned with the outflow which might be its 
host star. No X-rays are detected from the candidate ionizing star or from the jet. 

We thus find that the spatial distribution of the Chandra stars (Figure [8]) and their NIR 
counterparts (Figure [9]) show the following features: a compact cluster previously recognized 
as NGC 2237 which is ~ 2.5 x 1.3 pc in size, elongated along an east-west axis parallel to 
the direction of the NGC 2244 early O stars. Its shape might reflect an origin in a cometary 
cloud that was triggered to form stars and ablated by the O star ultraviolet photons and 
winds. A broader distribution of stars with a possible sparse subcluster is seen over a 5 x 5 pc 
region northwest of the main cluster, partly in the evacuated Rosette HII region and partly 
embedded in the optically dark pillars and cloud. The young star density is distinctly reduced 
1 — 3 pc south and west of the cluster; the western deficit can be attributed to absorption from 
the residual molecular cloud material. This NGC 2244 satellite cluster contains about ten 
intermediate-mass stars, consistent with its estimated population from the 2MASS study. 
There are hints of mass segregation (the intermediate mass stars lie primarily in the HII 
region) and age segregation (the K-hand excess stars like primarily along the cloud edge), but 
a clear age sequence is not seen. We note that richer substructures in the stellar distributions 



of m assive y oung clusters hav e been seen in our X-ray studies of N GC 6357 (IWang et al. 



20071 ). M 17 flBroos et al.ll2007f ). NGC 2244 (Paper I), and NGC 6334 flFeigelson et al.ll2009h . 



4.3. K-band and X-ray Luminosity Functions 



Following |lJ (120051 ) and IWang et al.l (120081 ) we create the apparent i^^-band luminosity 
function (KLF) for NGC 2237 with the spatially complete 2MASS data, subtracting the KLF 
of a control field representative of the background distribution of NGC 2237 and reddened 
with A y = 2 mag. Our KLF, s hown in Figure [T2| is consistent with the one derived for NGC 
2237 in Bonatto fc Bica J2009h . which turns over beyond the completeness limit of 2MASS 
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[Kg = 14.3 mag; ISkrutskie et al.ll2006[ ). It is well fit with a power-law distribution in the 
range 12 < Kg < 14.3 mag, giving a slope of ~0.6. Within uncertainti es, this is marginall y 
consistent with model predictions of ~0.4 for clusters of similar ages (iLada fc Ladalll995l ). 
A sig nificant excess of stars is seen as a "bump" at Kg = 11.5 in the KLF (see also Figure 
12 in iBonatto fc Bical 120091 ). indicating the presence of more intermediate- mass stars than 
expected from a simple extension of the KLF slope. This perhaps reflects the contribution 
from a secondary cluster in the field or foreground/background stars. For comparison, the 
KLF for the young, embedded cluster RMC XC (Paper II) is also derived and shown in 
Figure [T2l No excess of i^-bright stars is see n in the KLF and the slope is 0.39, in agreement 
with the KLF slopes in iLada &: Ladal (Il995l ). 



We further examined the 13 f^-bright stars to identify their nature. The contamination 
of bright stars from NGC 2244, the central cluster in the Rosette Nebula, is too small to 
account for the if-bright stars found. Radial profiles of NGC 2244 (Figure 2 in Li 2005; see 
also Paper I) show that the stellar surface density of NGC 2244 is approximately 0.2 star 
per square arcmin at a distance r ~ 15' from the cluster center and has already decreased 
to zero at r ~ 19'. At most two NGC 2244 stars may be present in the NGC 2237 region, 
which would not necessarily be i^-bright stars even if they are detected. 

Three of these i^-bright stars are not detected as X-ray sources, so they are most likely 
nearby foreground stars along the line of sight towards NGC 2237. For the remaining 10 
X-ray-detected K-band bright source^, their median photon energy ranges between 1 and 
1.5 keV. None shows very soft emission (median energy <^ 1 keV) that is suggestive of fore- 
ground stars with low absorption. Using B and V magnitudes from the NOMAD catalog 
( IZacharias et al.l 120051 ) . we derived the optical CMD for these stars, show n in Figure fT3l 
Adopt ing the tabulated spectral types and photometric properties of stars in ISchmidt-Kaler 
(119821 ). the zero-age main sequence (ZAMS) isochrone is shown, assuming a minimum red- 
dening of E{B — V)=0.5 and a distance of c? = 1.4 kpc for NGC 2237 members. Four 
stars (#54, 59, 83, 138) appear to be likely foreground FGK stars; the ZAMS location at 
d = 1 kpc with no reddening is shown for reference. The rest are consistent with being 
intermediate-mass stars in NGC 2237. 



Similarly we derived the observed hard band XLF (see lWang et al.ll2008l ) for NGC 2237, 
which is more robust against the inhomogeneous obscuration. Luminosities were taken from 
the spectral fits for the brighter X-ray sources, and estimated with a simple absorbed thermal 
plasma model for the faintest sources with few counts. We adopt an Nh = 4 x 10^^ cm~^ 
(derived from Nu/Ay = 1.6 x 10^^ cm~^ mag~^; Vuong et al. 2003, and a typical measured 



^These sources are ACIS# 44, 53, 54, 59, 71, 83, 128, 129, 138, and 141. 



- 14 - 



Ay ~ 2.5 ma g from Figure HI) , and a kT = 2 keV based on the median photon energy of the 
faint sources ( iFeigelson et aLll2005l ). The XLF with a farther distance d = 3.9 kpc was also 
generated similarly, except that a higher Nh = 5.5 x 10^^ cm~^ is adopted (Bonatto & Bica 
2009). 

Figure [H] shows the resulting XLF. The uncertainties for the hard band luminosity Lh 
derived from spectral fits range between 0.1 and 0.3 dex. To include both the y/N counting 
error and the Lx-dependent measurement error, following Getman et al. (2006), XLF error 
bars shown in Figure [T^ cover the range containing 68% (lex equivalent) of the Monte Carlo 
simulated XLF distributions, where individual X-ray luminosities are randomly drawn from 
Gaussian distributions with mean equal to the measured log Lh and variances ~0.3 dex 
(observed AlogL/^). 

For comparison, the X LF for NGC 2244 and the XLF for ONC from the Chandra Orion 
Ultradeep Project (COUP iGetman et al.ll2005bl ) were included. Considering log Lh > 29.6 



erg s~^, where the detections are mostly complete, the slope of the NGC 2237 XLF matches 
the ONC XLF well, but a scaling factor ~0.2-0.3 is required to fit the ONC X-ray-detected 
population . Given that the ONC contains ~2000 members with mass above the stellar 
limit (e.g., iHillenbrandl 119971 ) . this implies a smaller cluster population ~400-600 for NGC 
The X-ray-inferred population agrees well with the NIR population of ~230 stars 



le.g. 

2237. 

above 0.8 Mq, 



which is ~40% of the total population assuming a standard IMF. In addition, 
comparison of the NGC 2237 KLF to the RMC XC cluster KLF (Figure ^ also supports 
that NGC 2237 is less rich than the RMC XC cluster, which has an estimated population of 
~800 stars (Paper II). We conclude that NGC 2237 contains ~400-600 members. 



4.4. Discussion on the Distance to NGC 2237 



The adopted distance to NGC 2237 in this w ork is (i ^ 1.4 kpc, consistent with the gener- 
ally ac cepted distance to the Rosette region (e.g.. lOgura fc IshidalE^ : iRoman-Zuniga Sz Ladal 

20083k 

A spectroscopic survey is needed to securely establish the distance to NGC 2237. 



In iBonatto fc Bical (120091 ). the decontaminated near- infrared CMD is not well constrained 
and allows alternative age and distance solutions. For assumed ages of 1 Myr and 5 Myr, 
they derived be s t-fit d istances of d = 3.9 ± 0.4 kpc and d = 1.7 ± 0.3 kpc, respectively. 
Bonatto fc Bical (120091 ) favor the further distance with an age 5 ± 4 Myr since it better ac- 
counts for the decontaminated CMD sequence. NIR photometry of the B5 star GSC00154- 
02384 agrees with the larger distance [d = 4.3 kpc) while photometry of the A2 star 
GSC00154-01659 agrees with the closer distance {d = 1.4 kpc). We riote th at the spec- 
tral types were estimated from optical uvby photometry (jPhilip fc Egretlll980l ) and may be 
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uncertain. The A2 and B5 star are also detected in our Chandra observation as ACIS #128 
and #138, respectively. Both have ~ 10 net counts, and the full-band absorption-corrected 
X-ray luminosities logL^ c ~ 30 — 31 (for d = 1.4 kpc) are acceptable for inter mediate-mass 



stars, as the X-ray emission often comes from their lower-mass companions (e.g.. lStelzer et al. 



20051). 



Getman et al.l (j2008bl ) proposed an X-ray method to determine distance to young star 
clusters, based on two interlocking features: X-ray luminosities of PMS stars depend on 
stellar mass, as seen in a wide variety of PMS populations independent of the detailed age 
distribution and star formation mode (Preibisch et al. 2005; Telleschi et al. 2007); and X- 
ray luminosities and masses scale differently with the assumed distance. We attempted this 
method to obtain an independent estimate of the distance to NGC 2237 using a reliable 
sample of bright and less obscured {Nh < 10^^ cm~^) X-ray sources with c derived in 
Table |3] and mass M derived from NIR photometry and Siess et al. (2000) tracks. However, 
comparison with the well calibrated relation L^ c = 1.991ogM + 30.43 erg s~^ for the Cepheus 
0B3b and Orion stars (Getman et al. 2006) does not provide a decisive distance for the 
NGC 2237 cluster: a 2 Myr old cluster at d = 1.4 kpc and a 5 Myr old cluster at d = 3.9 
kpc both follow the relation with comparable scatter. 

The XLFs discussed in § 14.31 are more effective and support a nearer distance. Using 
d = 3.9 kpc, the NGC 2237 XLF implies that it is a rich massive cluster ~ 1.5 — 2 times 
the population of the ONC and NGC 2244. Dozens of OB stars should be present for a 
standard IMF, and these are absent in the optical, NIR, and X-ray images. In addition, the 
X-ray-measured absorption columns are moderate, a few times 10^^ cm~^, consistent with 
Ay ~ 1 — 3 mag seen towards the central Rosette NGC 2244 cluster. For a distance d = 3.9 
kpc along the Galactic plane, a significantly higher Ay is expected. 

Therefore, we consider our choice of the near distance justified. Note that the larger 
distance to NGC 2237 is not ruled out by the rich population estimated from the XLF com- 
parison; an anomalous IMF could explain the lack of OB stars. Further optical spectroscopic 
surveys of the NGC 2237 cluster will be invaluable to confirm whether it is a background 
cluster behind the Rosette region at a distance of d = 3.9 kpc. 



5. The Global View of Star Formation Across the Complex 
5.1. The Relation of NGC 2237 to NGC 2244 



NGC 2237 was dubbed a twin cluster (or satellite cluster) to NGC 2244 byllJ fl2005[ ). 
as the stellar density distribution of 2MASS sources shows a distinctive core to the west 
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of the compact core of NGC 2244 (see Figure 1 in Li 2005). iRoman-Zufiiga et al.l ( l2008a( ) 
also identified the cluster, and pointed out t hat the region was h istorically recognized as 
NGC 2237. The region was not covered in the iPoulton et al.l (120081 ) Spitzer survey. Notably 
the NIR-excess disk fractio n is significantly lower than those in embedded clusters (see 
Roman-Zufiiga et al.ll2008al and Paper II), but similar to that observed in the NGC 2244 
cluster. 

The NGC 2237 satellite cluster is very similar in various properties - obscuration, rela- 
tive location and disk fraction to the NGC 2244 cluster - to the X-ray-discovered RMC XA 
cluster on the eastern side of the Rosette Nebula (Paper II). The separation between NGC 
2237 and the NGC 2244 cluster core is ~14' (~6 pc), almost identical to the distance that 
separates RMC XA from NGC 2244. It is plausible that pre-existing massive molecular 
clumps in the natal cloud that formed NGC 2244 also collapsed and formed stars in RMC 
XA and NGC 2237 simultaneously with the formation of NGC 2244. However, they could 
also be triggered to form shortly following the onset of the HII region around the NGC 2244 
O stars. Either way, these three clusters represent the early star formation episode in the 
Rosette complex. The stars in the three clusters are probably indistinguishable in age, and 
they share a similar low circumstellar disk fraction (§ 13. 2p . NGC 2244 clearly has the dom- 
inant OB population; NGC 2237 has many lower-mass stars but few OB stars, and RMC 
XA is less rich. All three clusters show substructure: a subcluster around the 05 star HD 
46150 in NGC 2244 (Paper I); XAl, XA2, and XA3 in RMC XA (Paper II); and elongated 
structures in NGC 2237 ( 



Further towards the west of the satellite cluster, where the molecular pillars protrude 
into the HII region, the ionization front from the OB st ars in NGC 2244 ma y have triggered 
some star formation. These "elephant trunk" globules ( ISchneps et al.lll980f ) and bright rims 
(see Figure [T] and optical images of the Eagle Nebula) are commonly seen at the peripheries 
of many nebulae and are interpreted as illustrations of the massive star feedback into the 
parent molecular cloud. The evidence for triggered star formation in photoevaporating pillars 
is tantalizing. The alignment of the Class II sources in NGC 2237 may be a hint of an 
ionization front that triggered collapse along the boundary of the HII region (§ 14. 2p . but 
the number of stars involved is s mall. No clear evidence of an age gradient (or evolutionary 
sequence) is found. iDent et al.l (|2009[ ) have studied the nature and mass of the small CO 
clumps in the Rosette including the elephant trunks/globulettes in the NGC 2237 region. 
They show that the globulettes are mostly low mass, which are being probably accelerated 
by the OB winds, and their possible age is shorter than 1 Myr. The removal of the rim gas 
proceeds quite rapidly, thus NGC 2237 must have emerged relatively fast from its original 
clump. 
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Only a handful of new stars are forming at the vertices of the elephant trunks in 
NGC 22 37, perhaps fewer tha n the number of protostars at the tip of molecular pillars 
in M16 (iHester fc Deschl 120051) . The masses of the clumps in the expanding ring range 
between ~3-100 Mq (jPent et al.ll2009[ ). Star formation efficiency of a few percent will suf- 
f ice to form the stars observe d in the mass range of NGC 2237. Dynamical instability 
(IGarcia-Segura fc Francdll996l) appears to be not effective in forming stars here. As shown 



m lWhalen fc NormanI ( l2008l ). expanding HII regions can either stimulate or lead to turbulent 
flows against the formation of later generations of stars, depending on the efficiency of cool- 
ing within the shocked material. The observed deficiency of stars forming at the ionization 
front here may imply that the cooling is weak and dynamical instability is not promoting 
formation of second-generation stars efficiently. 

Note however that the paucity of X-ray stars in the pillars could be due to low sensitivity 
of our observation in this region. The X-ray luminosi ty of protosta r s may also be intrinsically 
low in the evaporating gaseous globules (EGGs): iLinsky et al.l ( 120071 ) report X-ray non- 
detection of EGGs in M16 and conclude that either the EGGs do not contain protostars or 
the protostars have not yet become X-ray active. 

More X-ray sources are seen distributed ahead of the ionization front; these account 
for the density enhancement northwest of NGC 2237 seen in Figure El Most of these X-ray 
sources show hard X-ray emission, suggesting that they are embedded young stars. The 
presence of a likely HH object (§ 14. 2p supports the ongoing formation of stars in this region. 
These stars may have formed as a result of interaction between the molecular materials here 
and the Rosette Nebula. Spectroscopy and Spitzer mid-IR data can determine their ages 
and evolutionary stages. 



5.2. Large Scale Star Formation in the Rosette Complex 



Figure [T^ shows the ROSAT PSPC image with complete spatial coverage of the region 
revealing the brightest X-ray-emitting stars. For a high resolution comparison, we assemble 
our Chandra images (exposure c orrected) in a ~ 1 .25° x 0.25° mosaic, as shown in Figure [T5b: 
a similar mosaic was shown by iTownsley et al.l (120031 ). The clumps of stars are apparent, 
which gives an overall view of their relative distributions, although the detection sensitivity 
is not uniform. To better understand the star formation history of this complex, we need to 
examine the large scale cluster formation events as an integrated process and look beyond 
our X-ray coverage. 



To outline the relative locations of the clusters with respect to the molecular material 
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and the morpholo gy of the ISM, Fig ure [T6k shows the clusters superposed on the ^^CO 
emission map from lHeyer et al.l (|2006[ ). The NGC 2244 cluster has apparently expelled most 
of the molecular gas while creating the central cavity. The RMC XA cluster and NGC 2237 
are mostly unobscured without surrounding molecular material, and the stars have low NIR- 
excess disk fraction. They may have formed nearly simultaneously or immediately after the 
formation of NGC 2244. The estimated total population is not as rich as NGC 2244, on 
the order of 200-600 (~10-30% of the total stellar population in NGC 2244). However, their 
stellar masses were probably determined by the available molecular material in their parental 
clumps; no strong feedback that may disrupt the cloud and suppress the star formation is 
expected from the NGC 2244 OB stars yet as they were forming. 

Moving further radially outwards, 25'-30' (10 — 12 pc) from the central OB clus- 
ter, a number of other clusters are present: RMC XB from Paper 11 (PL2, cluster C in 
Poulton et al. 2008, or structure B in Li & Smith 2005) to the southeast, PL 1 (clus- 
ter A in Poulton et al. 2008, also st ructure B in Li & Smith 2005) to the south, cluster 
REFLIO ( iRoman-Zuniga et al.ll2008al ) to the north. We also noticed a significant overden- 
sity of stars associated with the molecular clump to the west of NGC 2237 in both the 
2MASS study (see Figure 5 in Li 2005) and the FLAMINGOS NIR s urvey (Figure 27 in 
Roman-Zufiiga et al.ll2008a] ). The core is clearly seen in the CO map in iHeyer et al.l (|2006[ ) 
that shows prominent, dense molecular clumps further westwards from NGC 2237 (approx- 
imate center R.A. =06^^29™, Dec.=-|-05°00'). This cloud core wa s previously identified a s 
one of the principal r egions in the RMC, designated as Core E ( iBlitz &: ThaddeusI Il980l ). 
Roman-Zufiiga et al.l (l2008aj) find a collection of highly reddened NIR stars without J-band 
detections here. Neither Chandra nor Spitzer satellites have covered this molecular core. 
Nevertheless, we speculate that it could well be another region that harbors an undiscovered 
small embedded cluster. 



As seen in Figure [TBk. there is a prominent large shell-like structure in the CO map, with 
the above clusters located along the shell coinciding with the densest spots. Presence of such 
dense molecular material has been reported in other recent studies o f star formation tr iggered 
by the expansion of an HII region (e.g., iDeharveng et al.l l2009l : iKang et al.l l2009l ). The 



distribution of the clusters is also simplified in the cartoon presentation in Figure ITBb. These 
clu sters could be fo r med through the so-ca lled "coUect-and-coUapse" process (for review. 



see 



ElmegreenI Il998t IDeharveng et al.l |2005|) involving fragmentation of swept-up shells of 



interstellar gas by the expanding HII region of NGC 2244. When critical condensation is 
reached, the l ayers of neutral gas bec ome gravitationally unstable and eventually fragment 
to form stars ( IWhitworth et al.lll994j ). 



To elucidate this coUect-and-coUapse model, we adopt typical values in models for cloud 
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fragm entation induced by expanding HII regions and stellar wind bubbles (IWhitworth et al. 
1994J ). Assuming an isothermal sound speed = 0.2 km s~^ in the shocked gas, an emitting 
rate of hydrogen ionizing photons N^yc ~ 10^^ s~^ for the 04 and 05 stars, a wind luminosity 
~ 10^^ erg s~^, and a hydrogen nuclei number density ~ 1000 cm~^ in the interstellar 
medium, the timescale at which fragmentation occurs is ~ 1 — 2 Myr and fragmentation 
starts in the shell at radii of 6 — 10 pc with fragment masses around 10 — 23 M©. The ring 
of clusters around 10 - 12 pc (RMC XB = PL 2, PL 1, and REFL 10) correspond to this 
predicted ring of collapsing cloud fragments. Note the presence of a massive Class I protostar 
[Chandra source #89 = IRAS 06306+0437 in Paper II) in a sparse cluster within RMC XB. 

Several observations might confirm this coUect-and-collapse scenario. The model pre- 
dicts that active star formation, small clusters with high disk fractions, would be present in 
PL 1 and REFL 10 on the CO rim. The model also suggests that PL 2 formed more recently 



than ra ore distant clusters in the RMC. Consistent with this prediction, iRoman-Zuniga et al. 



(]2008al ) found that the IR excess fraction in PL 2 (33%±6%) is higher than PL 4 (24%±3%). 
The collect-and-coUapse process that formed PL2 may be enhanced from the expansion of 
cluster PL4 and PL5. 

We note that the four cloud fragments around the central NGC 2244 cluster form a 
bowl that is open towards the northeast, without significant molecular mate rial or clusters 



In th i s direction lies the large Monoceros Loo p supernova remnant (SNR; [Fountain et al. 



19791 : IVidal-Madjar et al.l Il982l : lOdegardI Il986[). and one m ay speculate that it may have 
some influence to the RMC. However, iTownsley et al.l (|2003[ ) reviewed the possibility of an 
encounter between the SNR and the molecular cloud, and there is no sign of interaction. 

Farther away from NGC 2244, the very young clusters PL4, PL5, REFL08 (together 
they are recognized as cluster RMC XC in Paper II, as structure C in Li & Smith 2005, 
and as cluster E in Poulton et al. 2008), PL6, PL7 (cluster G in Poulton et al. 2008), 
and REFL09 (cluster F in Poulton et al. 2008) are located along the densest mid-plane of 
the cloud. This was prev iously taken as evidence of sequentially triggered star formation 
( lElmegreen fc Ladalll977l ). where one cluste r provokes the format i on of t he next one. Based 
on the NIR studies of the embedded clusters, iRoman-Zuniga et al.l ( l2008a| ) propose a different 
idea of sequential formation. The IR-excess disk fraction among the clusters increases away 
from the nebula suggesting a rough sequence of ages, but 2 Myrs may not be enough time 
to produce a triggered sequence. The clusters PL7 and REFL09 are too close in age and 
distance to be triggered by clusters PL4 and PL5. Their formation is likely an outcome of 
the evo lution of the primordial mol ecular cloud, and not from the interacti on with the HII 
region ( IRoman-Zuniga et al.ll2008al ). In fact, recently IWilliams et al.l ( 120091 ) show that star 
formation in PL06 must have been active for about 2 Myrs and still continues. This reinforces 
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the view that cluster formation in the RMC may have proceeded without significant time 
difference, except that NGC 2244 has more OB stars with powerful winds and clears its 
surroundings more rapidly than the other low mass star clusters (e.g., NGC 2237, REFL 10) 
and deeply embedded small clusters (e.g., PL 06). 



6. Conclusions 

We present the first high spatial resolution X-ray image of the NGC 2237 cluster in the 
Rosette Nebula obtained via a single Chandra observation, which is part of our survey of the 
star formation activities in the Rosette complex. Prior to th is study, only 36 i^'-band excess 



cluster members were identified (IRoman-Zuiiiga et al.ll2008al ). The effectiveness of the X-ray 



survey technique is clearly demonstrated. Our main findings are: 

1. In this 20 ks observation, 168 X-ray point sources are detected with a limiting X-ray 
sensitivity of Lt ~ 10^° ergs s~^, nearly complete to the solar mass range. We associate 
134 of the 168 X-ray sources with cluster stars and infer about 24 other embedded sources, 
giving a total of ~ 160 stellar members identified here in NGC 2237. 

2. The locations of most ACTS sources in the NIR color-magnitude diagram are consis- 
tent with a small population of PMS low mass stars (M < 2Mq) in the NGC 2237 region 
with a visual extinction of Ay ~ 2 at 1.4 kpc, assuming a similar age to NGC 2244 (2 Myr 
old). Our X-ray sample provides the first probe of the low mass population in this satellite 
cluster to NGC 2244. We derive an overall ii'-excess disk frequency of 13% for stars with 
mass M > IMq using the X-ray- selected sample, consistent with the reported disk fraction 
from NIR study. 

3. We examine the X-ray properties of a number of B stars in the field. Several X-ray 
emitting stars are located at the peripheries of pillar objects. A previously unknown optical 
outflow feature is seen originating inside the optically dark region, which may be a newly 
identifled HH object; this supports the ongoing formation of stars in this region. 

4. The X-ray selected population provides a reliable probe of the cluster center and 
structure, which agree well with previous NIR studies. Similar to RMC XA, collapse of 
pre-existing massive molecular clumps accompanying the formation of the NGC 2244 cluster 
may have formed NGC 2237 as a satellite cluster. There is tantalizing evidence suggestive of 
a triggered formation process in the optical dark pillars by the NGC 2244 O stars, including 
discovery of young Class II sources aligned in an arc and embedded hard X-ray sources across 
the optical Photon Dominated Region. 
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5. We derive the KLF and the XLF for the NGC 2237 cluster. A significant excess of 
stars with apparent i^'^ ~ 11 — 12 is seen in the KLF, hkely associated with the presence of 
a few foreground stars. The slope of the NGC 2237 XLF matches the ONC XLF well, and 
the relative scaling implies a smaller cluster population ~400-600 for NGC 2237, which is 
consistent with the estimates from NIR studies. We discuss the uncertainty in the distance 
to NGC 2237 and prefer the near distance. 

6. The large-scale star formation in the Rosette Complex is reviewed in the context 
of the "coUect-and-collapse" scenario, given the locations of clusters and the morphology of 
the ISM. The formation of NGC 2244, RMC XA, and NGC 2237 represents the earliest star 
formation episode in the complex. The IRAS lOO/xm image and CO emission map show sites 
of currently active star formation-four potential clusters associated with the fragments of 
the large shell of swept-up neutral material from the expanding HII regi on. In the RMC 



a temporal sequence of star formation across t he complex is present (see iLi fc Smith! 12005 



Wang et al.l l2009l : iRoman-Zuniga et al.l l2008al ) . The cause of the sequence appears to be 



related to the molec ular cloud, perhaps primo rdial, not limited to the originally proposed 



triggering scenario in lElmegreen fc Ladal (119771 ) where one cluster provokes the formation of 
the next one. 
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Table 1. Primary Chandra Catalog: X-ray Source Properties 
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Note. — Table [T] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form 
and content. 

Note. — Column 1: X-ray catalog sequence number, sorted by RA. Column 2: lAU designation. Columns 3,4: Right ascension and declination for epoch 
J2000.0. Column 5: Estimated random component of position error, la, in arcseconds computed as (standard deviation of PSF inside extraction region)/ 
^(# of counts extracted). Column 6: Off-axis angle in arcminutes. Columns 7,8: Estimated net counts extracted in total energy band (0.5-8 keV) with 
Poisson-based standard errors. Column 9: Background counts extracted (total band). Column 10: Estimated net counts extracted in the hard energy band 
(2-8 keV). Column 11: Fraction of the PSF (at 1.497 keV) enclosed within the extraction region. Note that a reduced PSF fraction (significantly below 90%) 
may indicate that the source is in a crowded region. Column 12: Photometric significance computed as (net counts)/(upper error on net counts). Column 13: 
Log probability that extracted counts (total band) are solely from background. Some sources have Pb values above the 1% threshold that defines the catalog 
because local background estimates can rise during the final extraction iteration after sources are removed from the catalog. Column 14: Source anomalies: 
g = fractional time that source was on a detector (FRACEXPO from mkarf) is < 0.9 ; e = source on field edge; p = source piled up; s = source on readout 
streak. Column 15: Variability characterization based on the Kolmogorov-Smirnov statistic in the total band; a = no evidence for variability (0.05 < Pks)\ 
b = possibly variable (0.005 < Pks < 0.05); c = definitely variable {Pks < 0.005). No value is reported for sources with fewer than 4 counts or for sources in 
chip gaps or on field edges. Column 16: Effective exposure time i ks: approximate time the source would have to be observed on axis to obtain the reported 
number of counts. Column 17: Background-corrected median photon energy (total band). 



Table 2. Tentative Chandra Catalog: X-ray Source Properties 
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Note. — Table [2] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form 
and content. 



Note. — Column 1: X-ray catalog sequence number, sorted by RA. Column 2: lAU designation. Columns 3,4: Right ascension and declination for epoch 
J2000.0. Column 5: Estimated random component of position error, la, in arcseconds computed as (standard deviation of PSF inside extraction region)/ 
of counts extracted). Column 6: Off-axis angle in arcminutes. Columns 7,8: Estimated net counts extracted in total energy band (0.5-8 keV) with 
Poisson-based standard errors. Column 9: Background counts extracted (total band). Column 10: Estimated net counts extracted in the hard energy band 
(2-8 keV). Column 11: Fraction of the PSF (at 1.497 keV) enclosed within the extraction region. Note that a reduced PSF fraction (significantly below 90%) 
may indicate that the source is in a crowded region. Column 12: Photometric significance computed as (net counts)/ (upper error on net counts). Column 13: 
Log probability that extracted counts (total band) are solely from background. Some sources have Pb values above the 1% threshold that defines the catalog 
because local background estimates can rise during the final extraction iteration after sources are removed from the catalog. Column 14: Source anomalies: 
g = fractional time that source was on a detector (FRACEXPO from mkarf) is < 0.9 ; e = source on field edge; p = source piled up; s = source on readout 
streak. Column 15: Variability characterization based on the Kolmogorov-Smirnov statistic in the total band: a = no evidence for variability (0.05 < Pxs)\ 
b = possibly variable (0.005 < Pks < 0.05); c = definitely variable [Pks < 0.005). No value is reported for sources with fewer than 4 counts or for sources in 
chip gaps or on field edges. Column 16: Effective exposure time i ks: approximate time the source would have to be observed on axis to obtain the reported 
number of counts. Column 17: Background-corrected median photon energy (total band). 



Table 3. X-ray Spectroscopy 
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063027.33+045347.7 


11.5 


2.4 


21.9 






1.2 


53.6 (■ 


■■,+0.2) 


29.84 


29.7 


29.8 


30.1 


30.6 


12 


063028.13+050530.4 


20.7 


3.4 


21.2 


{■ 


■■,+0.6) 


2.1 


53.5 (- 


-0.3, +0.3) 


30.1 


30.1 


30.1 


30.4 


30.5 


14 


063028.66+050137.8 


32.6 


4.7 


21.9 






0.6 (-0.3, +0.7) 


54.3 (■ 


■■,+1.0) 


30.3 


29.6 


29.7 


30.4 


31.3 


16 


063028.95+050512.4 


10.3 


2.1 


20.4 






5.8 


53.2 




29.9 


30.1 


30.1 


30.3 


30.3 


17 


063029.28+050113.2 


9.8 


2.2 


22.1 






10.7 


53.3 




29.3 


30.3 


30.4 


30.3 


30.5 


19 


063029.51+050053.1 


11.7 


2.5 


22.1 


(■ 


■■,+0.5) 


1.1 


-5.8 




29.29 


30.5 


30.5 


30.5 


30.6 


20 


063030.14+045949.3 


14.1 


2.8 


23.0 


{- 


-0.3, +0.4) 


3.4 


54.3 (■ 


■■,+1.3) 




30.7 


31.1 


30.7 





^Cols. (l)-(4) reproduce the source identification, net counts, and photometric significance data from Table [Tl 

''Cols. (5) and (6) present the best-fit values for the extinction column density and plasma temperature parameters. In four cases (#19, 113, 162 and 163), a 
powerlaw spectral model gives a significantly better fit than a single-temperature thermal model. In these cases, column (6) gives the powerlaw index F and column 
(7) gives the normalization of the spectrum in italics. Col. (7) presents the emission measure derived from the model spectrum, assuming a distance of 1.4 kpc. 
Quantities marked with an asterisk (*) were frozen in the fit. Uncertainties represent 90% confidence intervals. More significant digits are used for uncertainties <0.1 
in order to avoid large rounding errors; for consistency, the same number of significant digits is used for both lower and upper uncertainties. Uncertainties are missing 
when XSPEC was unable to compute them or when their values were so large that the parameter is effectively unconstrained. Fits lacking uncertainties, fits with 
large uncertainties, and fits with frozen parameters should be viewed merely as splines to the data to obtain rough estimates of luminosities; the listed parameter 
values are unreliable. 

■^X-ray luminosities derived from the model spectrum are presented in cols. (8)-(12): (s) soft band (0.5-2 keV); (h) hard band (2—8 keV); (t) total band (0.5-8 
keV). Absorption-corrected luminosities are subscripted with a c. Cols. (8) and (12) are omitted when logNjj > 22.5 cm~^ since the soft band emission is essentially 
unmeasurable. 



Table 4. Stellar Counterparts 





X-ray Source 










Optical/Infrared Photometry 










Seq 


CXOU J 


USNO Bl.O 


B 


R 


I 


2MASS 


FLAMINGOS 


J 


H 


K 


PhCcFlg 


# 




ID 


(mag) 


(mag) 


(mag) 


ID 


ID 


(mag) 


(mag) 


(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


063007.57+050440.7 










06300735+0504423 


063007+050442 


14.12 


13.26 


12.89 


AAAOOO 


2 


063016.82+050452.2 


0950-0096116 


20.03 


17.53 


16.05 


06301680+0504518 


063016+050452 


14.52 


13.74 


13.42 


AAAOOO 


3 


063018.34+045953.9 






















4 


063020.40+050351.7 


0950-0096157 


20.41 




16.90 


06302046+0503528 


063020+050352 


15.04 


14.25 


13.89 


AAAOOO 


5 


063022.82+045538.2 












063022+045539 




18.36 


17.43 




6 


063023.75+050555.8 






















7 


063025.18+045941.0 


0949-0093548 


20.41 


18.04 


16.33 


06302519+0459412 


063025+045941 


14.80 


14.20 


13.86 


AAAOOO 


8 


063025.39+050540.5 


0950-0096215 


18.79 


16.75 


15.39 


06302550+0505402 


063025+050540 


13.89 


13.15 


12.81 


AAAOOO 


9 


063026.53+050240.9 


0950-0096223 




18.50 


17.32 


06302649+0502411 


063026+050241 


15.61 


14.74 


14.37 


AAAOOO 


10 


063026.78+045341.5 












063026+045344 




18.22 


17.45 





Note. — Table |4] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form 00 
and content. 

Note. — The errors of the BRI magnitudes are < 0.02 mag in general for stars V < 17 mag and become > 0.05 mag for fainter stars V > 19 mag (Berghofer 
fc Christian 2002). The me an photometric uncertainties in FLAMINGOS are 0.058 ± 0.012 , 0.064 ± 0.018 , and 0.056 ± 0.014 for J, H, and K, respectively 
dRoman-Zuhiga et al.ll2008ah . The typical errors are much smaller for brighter stars (e.g., trj = 0.02 mag for stars with J < 16 mag). For 2MASS photometric 
quality flags, see Cutri et al. (2003). 

Note. — Columns 1—2 reproduce the sequence number and source identification from Table [Tl and Table[2l For convenience, [OI81]=Ogura & Ishida (1981), 
[MJD95]=Massey, Johnson, & Degioia-Eastwood (1995), [BC02]=Bergh6fer &; Christian (2002) 



53- 


= [OI81] 14 = 


= [MJD95] 104; spectral type BIV; pmRA= 


11.0 mas/yr, pmDE=-2.8 mas/yr 


54 = 


= [OI81] 10 = 


= [MJD95] 108; spectral type B2V; pmRA= 


-2.3 mas/yr, pmDE=-11.9 mas/yr 


61_ 


=V539 Mon 


= [OI81] 13 =[MJD95] 110=MSX6C G206. 1821-02.3456; pmRA=2.8 mas/yr, pmDE=0.4 mas/yr 


71 _ 


= [OI81] 12 = 


= [MJD95] 102; pmRA=6.8 mas/yr, pmDE= 


=0.6 mas/yr 


128. 


= [OI81] 35 


= [MJD95] 471; spectral type A2:; pmRA= 


-0.8 mas/yr, pmDE=3.6 mas/yr 


138 


= [OI81] 36 


= [MJD95] 497; spectral type B5; pmRA=f 


3.5 mas/yr, pmDE=2.1 mas/yr 



= [MJD95] 498; pmRA=-3.0 mas/yr, pmDE=1.9 mas/yr 

= [BC02] 11; known X-ray source; log L2:(ROSAT/PSPC)=31.01 erg s'^; pmRA=0.6 mas/yr, pmDE=-12.6 mas/yr 
= [MJD95] 653; pmRA=-1.0 mas/yr, pmDE=-5.4 mas/yr 



I 

I— ' 

I 
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Fig. 1.— A 30' X 30' DSS2 i?-band image of the NGC 2237 region. The box outhnes the 
17' X 17' ACIS-I field of view centered on NGC 2237. The fields to the east outline the 
deeper exposures of the rich NGC 2244 cluster at the center of the Rosette Nebula presented 
in Paper I. 
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Fig. 2.— Chandra ACIS-I image of the NGC 2237 region (0.5-8 keV) overlaid with 168 
source extraction regions. The source extraction region is based on the local PSF, enclosing 
approximately 90% of the PSF at 1.49 keV. 
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Fig. 3. — Adaptively smoothed, exposure corrected ACIS-I image of the NGC 2237 region 
where red represents soft X-ray emission (0.5-2 keV) and blue represents hard X-ray emission 
(2-7 keV). Red sources are stars with emission dominated by soft counts. Sources near the 
field edge appear larger due to optical distortion of the Chandra mirrors. 
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Fig. 4. — NIR J — H vs. H — K color-color diagram for Chandra stars with high-quality 
JHK photometry. The green circles and red triangles represent sources with and without 
significant i^-band excess. The black solid and long-dash lines denote the loci of MS stars 
and giants, respectively, from Bessell & Brett (1988). The purple dash dotted line is the 
locus for classical T Tauri stars from Meyer et al. (1997), and the cyan solid line is the locus 
for Herbig Ae/Be stars from Lada & Adams (1992). The blue dashed lines represent the 
standard reddening vector with crosses marking every Ay = 5 mag. The weak-line T Tauri 
stars are located in the domain between the blue dashed lines indicated as the WTTS locus. 
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Fig. 5. — NIR J vs. J — H color-magnitude diagram using the same sample and symbols 
as Figure |H The purple solid line and dashed l ine are the 2 Myr isochrone and the 1 Myr 
isochrone for PMS stars from Siess et al. respectively. The dash dotted line marks 

the location of Zero Age Main Sequence (ZAMS) stars. The blue dashed lines represent the 
standard reddening vector with crosses marking every Ay = 5 mag and the corresponding 
stellar masses are marked. 
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Fig. 6. — (a) The X-ray spectrum and the spectral fit for the second brightest X-ray source 
in the field (#54). (b) The X-ray spectrum and the spectral fit for the brightest X-ray star 
(#149) in the field. See text for the best-fit model parameters. 
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Fig. 7. — The spatial distribution of X-ray sources that do not have matched IR counterparts 
(shown as blue crosses) and sources that have hard median photon energies {MedE > 2.0 
keV; shown as red circles). About ten hard X-ray sources without matched NIR counterparts 
are located inside the HII region cavity; these are likely AGNs. The rest, located inside the 
optical dark pillar region, are probably embedded young stars. 
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Fig. 8. — Map of the stellar surface density, in units of log stars per arcmin^, for NGC 2237 
sources with 6 or more net counts. It is smoothed with a 1.5' radius sampling kernel. 
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Fig. 9. — The spatial distribution of the X-ray sources superposed on the optical DSS 
image. The X-ray-selected stars are Class II (red diamonds) and Class III (green circles) 
based on their NIR colors. Some X-ray sources without NIR counterparts (blue plusses) are 
likely embedded st ars. The molecula r gas is outlined by the contours (red) from the ^^CO 
emission survey of iHeyer et al.l (120061 ). 
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Fig. 10. — Mass-stratified spatial distribution of NGC 2237 stars: M > 2Mq (red boxes), 
and M < 2Mq (blue crosses). Masses are based on identified spectral types or estimated 
from NIR photometry. 
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Fig. 11. — The KPNO Ha image of a region on the interface between the cloud and ionized 
nebula. Extraction regions of X-ray sources are shown in red. Two optically dark pillar 
objects are circled. Source #86 (only 4 net counts) appears highly obscured with a visual 
extinction of ~ 20 mag in the NIR diagrams, possibly a star embedded in or behind the 
dense pillar. A likely Herbig-Haro object is marked with a box. 



-43 - 



NGC 2237 KLF 



10.0 - 



1.0 



0.1 




10 12 
Ks (mag) 



100.0 - 



RMC-XC KLF 



10.0 



1 .0 - 



10 



14 



16 



Ks (mag) 



Fig. 12. — (a) The Kg-hand luminosity function of NGC 2237 after correction for reddening 
(assuming Ay = 2.5 mag) and background contamination. The KLF is found to match a 
power-law distribution with a slope of 0.6 (dotted line). The vertical dashed line indicates the 
2MASS Ks band completeness limit, (b) The KLF for cluster RMC XC (assuming Ay = 8 
mag; Roman-Ziiniga et al.2008a). The power-law slope here is 0.39. 
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Fig. 13. — The optical CMD for the 10 X-ray detected K-band bright sources, with ACIS 
source sequence numbers labeled. The solid line shows the zero-age main sequence (ZAMS) 
isochrone (Schmidt-Kaler 1982) assuming a reddening of E{B — V)=0.5 and a distance of 
d = 1.4 kpc, whereas the dashed line represent the ZAMS location for foreground stars at 1 
kpc with no reddening. The typical error bars for the photometry are shown. 



-45 - 



100.0 



10.0 - 



1 .0 



0.1 



2\ 



■ NGC 2237 (1.4 kpc) 
; NGC 2237 (3.9 kpc) 



29 



NGC 2244 
COUP ONC - 



L 



h 



30 



31 



32 



log Lx (ergs s" 



33 



Fig. 14. — The hard band (2-8 keV) X-ray luminosity function of NGC 2237 (black solid 
line). For comparison, the XLFs for NGC 2244 (red dotted line) and the COUP ONC (blue 
dashed line) are plotted. The XLF derived using a farther distance d = 3.9 kpc (green solid 
line) is also shown. 



Fig. 15. — (a) ROSAT PSPC image of the Rosette Complex with existing Chandra pointings 
outhned by the red boxes, (b) Chandra ACIS-I mosaic of six observations in the Rosette 
Complex. Red represents soft band (0.5-2 keV) X-ray emission, and green represents hard 
band (2-7 keV) X-ray emission. The Chandra image is rotated for better layout in the page; 
the actual orientation is the same as (a), where nort h is up and east is to t he left. A similar 
image without the NGC 2237 exposure is shown by iTownsley et al.l (120031 ). 




Fig. 16. — (a) A large scale view of th e Rosette Complex using the lOO/xm IRAS map shown 
in color. ^^C O emission contours fro m Hever et al. I J2nn6h a re overl aid. Previously identified 
clusters from IPhelps fc Ladal (119971 ) , iRoman-Zuiiiga et al.1 ( ]2008al ) , and IWang et al.l (120091 ) 
are labeled, (b) A cartoon presentation of (a) with the dark clouds and nearby supernova 
remnant outlined and the clusters marked with stars. 



